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We investigate the spin-induced polarization in the multiferroic compound Ba2CuGe207 using 
variarional and finite temperature mean field approaches. The compound is described by a spin-3/2 
Heisenberg model extended with easy plane anisotropy and Dzyaloshinskii-Moriya (DM) interaction. 
Applying magnetic field parallel to the [110] axis, three phases can be distinguished: (i) At high 
magnetic field we find a partially magnetized phase with spins parallel to the fields and uniform po- 
larization; (ii) Below a critical field the ground state is a twofold degenerate canted antiferromagnet, 
where the degeneracy can be lifted by a finite DM interaction; (iii) At zero field a (7(1) symme- 
try breaking phase takes place, exhibiting a Goldstone-mode. To reproduce the magnetization and 
polarization measurements reported in Murakawa et al. [Phys. Rev. Lett. 105, 137202 (2010)], 
we introduce an additional term in the Hamiltonian that couples the polarizations on neighboring 
tetrahedra. This results in the appearance of a canted ferrimagnetic phase for h < IT, characterized 
by a finite staggered polarization, as well as by a finite magnetization along the [1, 1,0] axis that 
leads to torque anomalies. 

PACS numbers: 75.85.+t 75.30.Gw 75.10.-b 



I. INTRODUCTION 

In the traditional sense multiferroic materials are 
characterised by the coexistence of ferroelectricity and 
ferromagnetismii"— In these "proper" multiferroics the 
magnetoelectric effect is quite small, since the simulta- 
neous presence of the two order is rather difhcult due to 
the fact that they break the inversion and time rever- 
sal symmetry in a different way: while the ferroelectric 
order breaks space inversion symmetry and is invariant 
under time reversal, the magnetic order behaves just in 
the opposite way. In addition, the concurrent appear- 
ance of these orderings does not necessarily mean that 
the magnetic and electric dipoles are strongly coupled to 
each-other. 

After almost fifty years of research in the field, the dis- 
covery of the giant magnetoelectric response in TbMnOg^ 
has launched a new concept: spin driven ferroelectricity. 
In the last decade - along with experimental realizations 
- many theoretical explanations have been proposed as 
the source of electric polarization in magnetically ordered 
materials. Electric polarization induced by cycloidal 
spin order was explained through spin chirality^ or in- 
verse Dzyaloshinskii-Moriya mechanism^ and was found 
experimentally in many materials, including TbMnOa;^ 
NiaVaOg,^ CuFeOa,^ MnW04,^ CoCr204,^° LiCuaOa" 
and CuO^^, just to mention a few. Exchange striction 
can also lead to electric polarization, as was shown in 
the case of the perovskitc i?Mn03 materials, with R be- 
ing a rare earth ion.^'^ The aforementioned mechanisms 
may induce polarization jointly, as predicted in the case 
of i?Mn205 materialsi^^-i^ In all these cases the mecha- 
nism involves two spins on a bond. 

If the spin is located on a site that breaks the inver- 
sion symmetry, the polarization can be induced by a sin- 



gle spin. The origin of this process is the spin depen- 
dent metal-ligand hybridization arising from spin-orbit 
coupling liiii^ Murakawa and collaborators proposed that 
this mechanism explains the induced ferroelectric polar- 
ization in Ba2CoGe207J^ 

As a result of strong easy-plane anisotropy, below 
T/v = 6.7 K the S = 3/2 moments order into a canted 
antiferromagnetic pattern that is confined in the Co- 
plane. ^'i This canted planar antiferromagnetic phase is 
a multiferroic phase in which magnetoelectric behavior 
has been observed. Ascribed to the symmetry properties 
of Ba2CoGe207 the sum over the vector spin chirality 
Si X Sj vanishes and the exchange interaction Si • Sj is 
uniform for all bonds. Hence the induced polarization 
cannot be explained by the concept of spin current or 
exchange striction. Spin dependent hybridization mech- 
anism was suggested as a solution, according to which the 
local polarization takes the form of P oc X]i=i(S ' 6i)^ei, 
where ei vectors point from the Co ions toward the sur- 
rounding four O ions.^^ This model intrinsically possesses 
the symmetries of the lattice and gives the same result 
for the magnetization dependence of the polarization vec- 
tor as a thorough symmetry analysis (e.g. in Ref. plj) 
would. It recovers the sinusoidal response of electric po- 
larization to the rotating magnetic field and describes the 
nature of induced polarization in magnetic field qualita- 
tively well. IS 

In this paper we investigate the magnetic field depen- 
dence of induced electric polarization in Ba2CoGe207. 
The paper is structured as follows: we start with a short 
symmetry analysis in Sec. |TT]to deduce the Hamiltonian 
and the spin dependent expression of the polarization. 
In Sec. mil we study the field dependence of polarization 
at zero temperature variationally, using a site factorized 
wave function. Two cases are considered, when the ex- 
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ternal field is applied along the [110] and [100] axes. We 
investigate the efi'ect of DM interaction and an additional 
antiferro polarization term in detail. In Sec. IIVI we per- 
form finite temperature mean field approximation to cal- 
culate the polarization and magnetization for both field 
settings. We qualitatively compare our ft,||[110] results to 
the experimental findings of Ref. [Toj and give predictions 
for the induced polarization with regard the /i||[100] field 
setting. 



II. SYMMETRY CONSIDERATIONS 

A schematic figure of Ba2CoGe207 is shown in Fig. [T] 
The Ba2CoGe207 has a tetragonal non-centrosymmetric 
structure characterized by the space group P42im. The 
neighboring C0O4 tetrahedra are tilted from the [110] 
crystallographic direction by the angles ±k, and due to 
the differently oriented tetrahedral environments the unit 
cell contains two Co^+ ions with 5 = 3/2. 

The point group of the lattice is isomorphic to 
the group I?2d that consists of the groups ^4 — 
{E, 5*4, C2, 5*4}, with the rotation axis located on the Co 
sites, and of C2-U = {i?, C2, crp-j^gj, cr[iio]} in the center of 
the four Co sites, as shown in Fig. [Ija). The symmetry 
properties of the lattice determine the terms that can be 
included in the Hamiltonian, i.e. the terms that trans- 
form as the fully symmetric irreducible representation 
Ai of I?2(i- When classifying the operators we assume 
that the magnetic order is a two-sublattice order, i.e. 
the magnetic unit cell corresponds to that of the lattice. 
Our Hamiltonian has the form of 
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-5h^s. + i?, {sfsy-sfsj), (1) 

where denotes pairs of nearest neighbor sites. The 
orientation of the DM vector (see Fig.[TJb)) is fixed by the 
point group of the unit cell. The lattice symmetries allow 
for other terms, too, such as the g-tensor anisotropy: 
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with hx — /i[iio] and hy = h^-^^^y Although, we found 
that it has no significant effect on the field dependence 
of the polarization, therefore we did not include it in our 
model. 

As mentioned before, the polarization P and the spin 
behave in different ways under inversion (T) and time 
reversal (T) operations: XP — and TP = P, while 
IS — S and TS — —S. For there is no inversion center in 
Ba2CoGe207, we only have to account for T, and it fol- 
lows that the polarization can be coupled linearly to op- 
erators that are even-order in spin and consequently are 
invariant under time reversal. With respect to the lat- 
tice symmetry, the z-component of polarization vector ~ 




FIG. 1: A schematic figure of Ba2CoGe207. (a) the point 
group of this material can be constructed from the groups S4, 
and C2v acting on the Co sites and in the middle of four sites, 
respectively. The mirror planes are labeled by their normal 
vector, (b) due to the tilted tetrahedral environment the unit 
cell contains two Co^^ ions that we denote by A and B. The 
symmetry allowed DM vectors - with respect to the order of 
involved sites (black arrows) - are indicated by purple arrows. 
At the r point the in-plane components of DM vector cancel 
each-other thus in our case it is sufficient to take only the 
out-of-plane components into account. 



or in fact any polar vector - transforms as the irreducible 
representation B2 of 2?2d, while the in-plane components 
(p^,p^) belong to the two dimensional representation E. 
Classifying the second-order spin operators according to 
lattice symmetry we find that, e.g. (S*^)^ — (S^)'^ and 



j j transform as the irreducible representa- 
tion B2 as well, therefore they can be coupled to the 
z-component of polarization. Similar logic leads to the 
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form of the operators and P^: 



parallel to the applied field h[ioo]- % = 1/V2{SJ ~ S'J) 
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{Syf - {S^f) - sin 2k, [s^S" + S'j'S'j)3) 



where j belongs to either sublattice A, or B. The dif- 
ferent orientation of the tetrahedra is accounted for by 
choosing Kj^A — and Kj^b = — Note that the oper- 
ators defined in ([3]) are actually quadrupole operators. 



III. ZERO TEMPERATURE VARIATIONAL 
APPROACH 

Based on the neutron scattering measurements of 
Ref. 20] that suggested a two-sublattice canted antifer- 
romagnet as the ground state, we chose a site factorized 
variational wave function: 
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where the |3/2), |l/2), ...denote the \S^) of the spin-3/2 
on the A and B sites. There are 6 independent variational 
parameters for \iJja) and another 6 for \^b) to be deter- 
mined by minimizing the ground state energy 



E = 
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For simplicity we assume that there is no exchange 
anisotropy: J = - as in Ref. [22|. We study the effect 
of the DM interaction in detail. 



A. Mapping between the cases /i||[110] and /i||[100] 

In the following sections we compare two cases: when 
the field is applied along the axes [110] and [100]. 
Whereas regarding the magnetization it is natural to con- 
sider m[iio] and its rotated counterpart m^ioo] as the rel- 
evant magnetic order parameters, it is not that easily 
perceived which component of the polarization will char- 
acterize the electric order. Forgetting about the polariza- 
tion, the two kinds of field settings can be mapped to each 
other by a 7r/4-rotation of the local basis. We introduce 

the new spin operators so that is the spin component 



1/V2{S}I + Sf). In the case /i[[[l,0,0], the 



Hamiltonian in the new basis has exactly the same form 
as HI) when the field is parallel to [1,1,0]. The polariza- 
tion operators, on the other hand, change. Considering 
the z-componcnt of P, in the new basis it will have the 
form of 
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with k = K + 7r/4. In order to an easier understanding 
of the relation between induced polarization in the two 
basis we sketched the effect of local basis transformation 
in Fig. [21 As it was pointed out in Ref. [l^, the induced 
polarization has extremum when the spin is aligned with 
either the upper or the lower edge of the surrounding 
tetrahedra. (The lower and upper edges are those paral- 
lel to the Co-plane.) If is maximal for a given spin 
orientation, the 7r/2-rotation of the spin will cause the 
polarization to change sign. When the spin lays halfway 
between these two positions, the polarization is zero. Ap- 
plying the field in the direction [110], the spins on the two 
sublattices are tilted (by the same angle) from the same 
symmetry plane of the tetrahedra, therefore the induced 
polarization on A and B has the same sign (and magni- 
tude). However a 7r/4 rotation of the local basis rotates 
the tetrahedral environment in a way that the spins on 
A and B are now tilted from the different edges, and as 
a result the induced polarizations have different sign on 
the two sublattices. (see Fig. ^ 



(a) 



[110] 





FIG. 2: Relation between the cases /i||[110] and /i||[100]. In 
our notation Pj takes maximum value for a spin aligned paral- 
lel to the lower edge (dashed line) and minimum value along 
the upper (solid) edge. (a)Schematic spin configuration for 
/i||[110] and Dz < 0. The spin induced polarization is the 
same on the two sublattices for the spins are canted by the 
same angle from the same (lower) edge of the tetrahedra. 
(b) Spin configuration for /i||[100] and Dz < after rotat- 
ing the basis by the angle 7r/4. Note that due to the rota- 
tion of the tetrahedral environment the induced polarization 
changes. The spins of the two subblatices are now canted from 
the different edges, consequently the induced polarization on 
them have different sign. 



According to this argument, when the field is parallel 
to [110] the relevant order parameter is the total polar- 
ization P^ = P^ + Pbi whereas for ft,||[100] the staggered 




- A 




component P| 



P^ — P\ will give a finite expecta- 



tion value. We shall note, that the difference between 
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the two cases is related to the lattice symmetries as well. 
While along the [110] direction the symmetry operation 
that connects site A with site B is a crjuo] mirror plane, 
along the [100] direction the sites can be transformed into 
each-other by a 2i screw axis, (see Fig. [I](a)) 

B. = 

Let us begin our investigations with the case = 
0. In zero magnetic field the variational ground state of 
([1]) is a planar antiferromagnet (superfluid) that breaks 
rotational symmetry ?7(1), and in which the staggered 
in-plane magnetic order parameter 0[/(i) = ^l^;^ — 

has finite expectation value, where S^^ — {S^ ,S^)?^ The 
planarity is the consequence of the easy-plane single ion 
anisotropy and the U{1) symmetry breaking is related to 
the fact that in the absence of an in-plane magnetic field 
the Hamiltonian commutes with total . The ground 
state wave function of site A can be expressed as 

1*^) = e-'^P^^^^f) (9) 
with a single variational parameter rjA: 

I^SF^ ^ If) + 1-1) +^iVA{\\) + \-\)) ^^^^ 

and a similar expression stands for l^'s) with parameters 
ifB and r/B- and ipB determine the angles of the spins 
with respect to the [110] direction. Furthermore we find 
that r]A = TjB = ?7 > 1 for the whole magnetization 
process. In zero field = and ipB = ^ — t^-, i-e. the 
spins are antiparallel;23i The t/(l) symmetry breaking is 
manifested in the fact that ip can take arbitrary value. 

In a finite magnetic field ip is not arbitrary any longer 
and the canted spins turn so that the total magnetiza- 
tion points into the direction of the field. This two sub- 
lattice canted order is twofold degenerate because, in the 
absence of the DM term, the Hamiltonian is invariant 
under the exchange of sites A and B. With increasing 
field the angle Sp = ipA — Pb between the spins of the 
sublattices decreases from Sp — tt to Sp — 0, while the 
length of the spins is unchanged (rj is constant). 

At a critical value (he) the two spins become aligned: 
PA ^ 'Pb = for the field parallel to [110] and pA = 
Pb = — 7r/4 for /i||[100]. The two sublattice order van- 
ishes and a uniform phase appears. The critical field, 
however, is not equal to the saturation field, the spins 
are not fully magnetized yet. With the further increase 
of the field the magnetization increases slowly as the spins 
reach their full length (i.e. ij decreases to = 1). 

The polarization and magnetization for Dz — in mag- 
netic fields along the axes [110] and [100] are plotted with 
black solid lines in Fig. ^c) and|4l[c), respectively. The 
spin configurations of the twofold degenerate canted or- 
der are shown in panel (a) and (b) of Figs. |3] and H) In 
both cases we indicated the canting angle of the spins 



with respect to the applied field. Conveniently, in the 
case of /i 1 1 [110] this canting angle is pA and pb ^ follow- 
ing from the definition of the ground state However, 
when /i]][100] we need to introduce the canting angles pA 
and Pb, to measure the angle from the [100]-axis. Actu- 
ally Pa and pb play the same role in the rotated basis 

introduced in Sec. IIII Al (when is the spin component 
along [100]) as pA and pb '^^ the original basis. Fur- 
thermore one could express the ground state ^ in the 

rotated basis, using pa/b = Va/b + ^/4. 



C. Finite 

When the DM interaction is finite, the angle between 
the spins of the sublattices A and B varies from tt, 
{pA — Pb 7^ TJ"), and the spins become canted even in 
zero magnetic field, with a canting angle depending on 
the magnitude of Dz- Nonetheless, this canted antifer- 
romagnetic state breaks the U{1) symmetry too, for Dz 
controls only the direction of the spins on different sub- 
lattices relative to each other, and remains a good 
quantum number (['^,5'^] = 0) as long as there is no 
field applied in the plane. A finite magnetic field lifts 
the f7(l)-degeneracy and we enter into the canted an- 
tiferromagnetic phase, where the canting angle depends 
on Dz and the field. The finite DM coupling - being 
sensitive to the exchange of A and B - lifts the twofold 
degeneracy of the canted phase and depending on its sign 
either of these two states is preferred. (A sign difference 
in D • {Sa X Sb) is equivalent to the exchange of sites A 
and B). Similarly to the case Dz — 0, the variational pa- 
rameters are rjA = TjB = rj, and below the critical field rj - 
consequently the spin length - is constant, as opposed to 
h > he, when the spins are partially magnetized and their 
length increases with the field. Comparing our results to 
the measurement of Ref. [l^ in the case of /i]][110], we 
find that Dz < is the appropriate choice. 

A schematic figure of the ground states in the canted 
phase is shown for /i]][110] in Figs.[3l^a) when Dz > and 
[31[b) when Dz < 0, and correspondingly in Figs. |ll[a) and 
(b) for h\\ [100] . The magnetization and the spin order in- 
duced polarization is plotted in Figs.ISjc) and|4l^c). Based 
on the argument of Sec. IIII Al when the field is pointing in 
the direction [110] the total polarization is finite, while for 
/i]][100] the staggered polarization gives non-zero expec- 
tation value. We remark that with the choice of k = 7r/8, 
the total polarization for /ij] [[1, 1, 0] and the staggered po- 
larization for /ij] [1,0,0] would match perfectly. A finite 
Dz interaction chooses different branches of the electric 
order parameter {Pz and Pf) for the two settings of mag- 
netic field, e.g. Dz < favors the higher branch of the 
total polarization for ft,]][l, 1,0] and the lower branch of 
the staggered polarization in the case of /ij] [1,0,0]. 
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FIG. 3: (a),(b) Schematic figure of ground states in canted 
antiferromagnetic phase when /i||[110]. For = (a) and 
(b) are degenerate. A finite DM coupling hfts the twofold de- 
generacy, the ground state configuration for > is shown 
in (a) while (b) is selected when D^, < 0. (c) The polarization 
Pz, the magnetization m^nQj and the staggered magnetization 
"^[Tio] ^® shown. Black line indicates the twofold degeneracy 
of the Dz = case below the partially magnetized phase. The 
colored lines correspond to the > and < cases in 
accordance with the coloring of spin states in (a) and (b). 




FIG. 4: (a),(b) the canted antiferromagnetic ground state for 
h\ I [100] . Similarly to the case oi h\\ [110] (a) and (b) are degen- 
erate when Dz = 0. Whereas a finite lifts this degeneracy, 
and Dz > selects the canted state (a) while Dz < prefers 
the configuration of (b). (c) When /i,||[100] the uniform po- 
larization vanishes (Pz — P^ + Pz = 0) and the staggered 
polarization Pf" = Pz — P^ gives a finite value. The coloring 
of the cases Dz > and Dz < follows that of the spin states 
used in (a) and (b). 



D. Effect of antiferro polarization term 

For higher fields the polarization curve is satisfyingly 
similar to the measurements, however, the low field be- 
havior cannot be explained using only the Hamiltonian 
([1]). In order to reproduce the sharp decrease of polar- 
ization below ft, 1 T we add an anitiferro [Kz > 0) 
polarization term to our model: 

Hpoi^KzY^PtPj (11) 

This is a kind of anisotropic biquadratic term that is al- 
lowed by the symmetry. Since [HpoijS^] ^ 0, the U{1) 
symmetry is lost with this term in the Hamiltonian. In- 
stead, a fourfold degenerate ground state appears (ac- 
tually this corresponds to the ground state discussed in 
Ref. i24i. where the order parameters in Ba2CoGe207 were 
investigated from purely the aspect of symmetries). We 
start our investigation with — again. The spin di- 
rection is determined by the minimization of (jlip and ([Ij . 
As it was discussed in Ref. [3 and Sec. HTTM the ori- 
entation of the spins relative to the surrounding tetrahe- 
dron determines the induced polarization, (see Fig. 3(a) 
in Ref. [HI). When the spin is pointing along the lower 
(or upper) edge of the tetrahedron, the induced is 
maximal (or minimal). Therefore a term P^Pb f^^vors a 
spin configuration in which the spins on sites A and B 
are parallel to different edges. (Note that this is not an 



orthogonal spin configuration for the tetrahedra of the 
different sublattices are rotated by 2k ss 48° compared 
to each other.) The polarization term competes with the 
antiferromagnetic exchange interaction and the resulting 
ground state is the canted state shown in Fig. [5{a). A 
finite magnetization points along the [100] or [010] axes 
in these cases. 

When the DM interaction is finite this spin configura- 
tion is favored hy Dz > that only changes the canting 
angle. However, a Dz < introduces frustration, as it 
prefers a reversed spin orientation: if the angle between 
the spins is ipA — for Dz > 0, it is —tpA + Vb for 
Dz < (see Fig. [51(b)). Loosely speaking this corre- 
sponds to the exchange of sites A and B, however with a 
smaller canting angle due to the competition of Kz and 
Dz. 

Selecting Dz < — that recovers the experimental re- 
sults — we calculated the induced polarization for both 
orientations of the magnetic field. The variational and 
the finite temperature mean field results are shown to- 
gether in Fig. [6] and [7] for ft,||[110] and in Fig. |8]when 
/i||[100]. In the former case we observed a new phase 
below hc2 ~ 1 T. This intervening phase is twofold de- 
generate. Following from the different spin length on 
the sublattices A and B a suitable order parameter is 
|Sa| — |S_b|, therefore we refer to this phase as a canted 
ferrimagnetic phase. The ground state can still be writ- 
ten as in (fTO|l. but here 'qA Vb- The angle between 
the spins on sublattices A and B is determined by the 
collective effect of Kz, Dz and /i[iio]. A schematic spin 
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a h=o 



C) 0<h<h„ 




[100] 




(d) 0<h 




FIG. 5; Ground state spin configurations when the antiferro 
polarization is present in the Hamiltonian. (a) and (b) the 
fourfold degenerate ground state in zero field for difTerent 
signs of . A spin configuration - with smaller canting an- 
gle - shown in (a) would correspond to the case Dz = 0, too 
(c) when /i||[110] a canted ferrimagnetic phase emerges be- 
low hc2 ~ 1 T due to the polarization term, (d) for /i||[100], 
since > 0, the canted ferrimagnetic phase is missing, and 
at finite field the ground state is the canted antiferromagnet 
similar to the case = 0. A ferromagnetic polarization term 
would have the opposite effect. 



configuration of tliis phase is shown in Fig. [SJ^c). When 
h > hc2 ~ 1 T we enter the canted antiferromagnetic 
phase that was characteristic in the case Kz = as well 
(see Fig. [3] (a)). When the field exceeds he ~ 13 T, the 
partially magnetized uniform phase emerges. In the ro- 
tated field setting {h\\ [100]) the intervening phase is miss- 
ing, and only three phases are observed. The finite field 
lifts the fourfold degeneracy of the h = Q ground state 
and we enter the non-degenerate two sublattice canted 
antiferromagnetic order shown in Figs. [Sl^d) or|3Ub). At 
the critical field the canted antiferromagnet is replaced 
by the uniform phase. We shall note that a ferro po- 
larization coupling [Kz < 0) would reverse the situation, 
and it would cause the emergence of the canted ferrimag- 
netic phase when the field is ft.||[100], while for /i||[110] 
this phase would not be present. 



IV. FINITE TEMPERATURE MEAN FIELD 
THEORY AND COMPARISON WITH THE 
EXPERIMENT 

Assuming site factorized solution for the finite tem- 
perature mean field calculations as well, we can write 
the Hamiltonian as "Hmf = T~iA + T~Lb, where 



-Ha = 4JS^(Sj 



ADz [SliSD - SliSD) 



and we obtain %b by exchanging the sublattices A and B 
and the sign of Dz ■ 'Ha and % b acts only on sublattice A 
and B, respectively and the factor 4 corresponds to the 



coordination number. The (S^) and 
thermodynamical average given as 



fP|) denote the 



(Sb) 



Tr(Si3e-'^^«(<^'^>'<-f'-J>') 
rpj.g-/3-HB((s^>,(Pl>) 

Tr(P|e^^'"«(<^'^>'<-f'-^>)) 
-^^g-mB{{SA)APl)) 



(13) 



(12) 



with (3 ~ 1/T. A similar expressions stand for (S^) 
and {PX) with the interchange of A and B. The self- 
consistent solutions of this set of equations, that consti- 
tute the mean-field approximation, is obtained by itera- 
tion. 

In our finite T calculations we took a realistic param- 
eter setting: J = = 1.885 K, A = 15.08 K according 
to Ref. [22], g = 2.2 and the values Dz and Kz were 
set to —0.02 K and 0.01 K respectively. The mean field 
result for T = 2 K is essentially undistinguishable from 
the T=0 variational calculation. For the field applied 
parallel to [110], we find that below hc2 ~ 1 T the polar- 
ization drops to zero and in this region there is a finite 
expectation value of m^jiigi and P|* corresponding to the 
zero temperature canted ferrimagnetic phase. The fact 
that these functions can take positive and negative values 
reflects the twofold degeneracy of the ground state. At 
higher field we enter the canted antiferromagnetic region 
and at there is a continuous phase transition to the 
partially magnetized uniform phase where the spins are 
aligned with the external field. With further increase of 
/i[iio] the spins grow to reach their full length. The finite 
temperature results for magnetic field /i||[l, 1, 0] are sum- 
marized in Fig. [6l Taking strictly the Hamiltonian ([1]) 
extended with the polarization term ()11|) we can recover 
the main characteristics of the experimental findings in 
Ref. fl^ . The polarization drops sharply below 1 T and it 
changes drastically with increasing temperature, whereas 
the magnetization curve is almost unchanged. Torque 
measurements would be conclusive to our studies, and 
could reveal if there is an orthogonal component of mag- 
netization at low fields {h < hc2)- Nonetheless, there 
are yet some properties to account for. The polariza- 
tion curve in the experiments starts from negative value 
at zero field, while in our model — due to the antiferro 
polarization term — Pz — when h = 0. Furthermore 
around hc2 the shape of the polarization is softer, in con- 
trast with our findings, that exhibit an edge in Pz when 
the canted ferrimagnet transforms into the canted anti- 
ferromagnetic phase. A possible solution to these anoma- 
lies can be the presence of a small orthogonal field. In 
Fig. [7|we plotted our results for the same setting as in 
Fig. [6l but with a small applied field in [TlO]. Naturally 
in this case the twofold degeneracy of the canted ferri- 
magnetic phase is lifted as the perpendicular field selects 
between ztm^yio]- 

In Fig. [5] we show the results for [1,0,0]. As men- 
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tioned before the relevant quantity in this case is the stag- 
gered polarization P|' . In good agreement with our zero 
temperature findings, the canted ferrimagnetic phase is 
absent, and below he the canted antiferromagnetic phase 
takes place. Similarly to /i||[110] at higher fields we reach 
the partially magnetized uniform phase. 




h[iio]n h[iio][T] 




5 10 15 20 



h[iio][T] 



(b) (c) 




1 2 3 1 2 

h[iio]n h[iio][T] 



FIG. 7: As a possible solution to the differences in our results 
and the measurements, we applied a small magnetic field per- 
pendicular to [110], e.g. /lelllTlO] and found that the char- 
acteristics of the experimental results are qualitatively well 
reproduced.(J = = f .885 K, A = 15.08 K, = -0.02 K, 
= 0.01 K and g = 2.2) 



FIG. 6: (a) behavior of magnetization and polarization in ex- 
ternal field along [110] at various temperature values, (b), (c) 
for fields /i ^ 1 T we observe the canted ferrimagnetic phase 
that is the consequence of the antiferro polarization term (|ll|l . 
In this phase the ground state is twofold degenerate as seen 
in the nature of Pf and rrifj^g^ .{J = = 1.885 K, A = 15.08 
K, = -0.02 K, = 0.01 K and g = 2.2) 



V. CONCLUSION 

We investigated the field and temperature dependence 
of the induced polarization in the multiferroic compound 
Ba2CoGe207. A detailed analysis has been given on the 
effect of DM interaction and an additional polarization- 
polarization term for two orientations of magnetic fields: 
/i II [110] and /i 1 1 [100]. We found that in the former case 
an antiferro polarization coupling leads to the emergence 
of a canted ferrimagnetic phase when < h < hc2 ~ 1 
T. In this region the polarization decreases sharply re- 
producing qualitatively well the experimental findings in 
Ref. [19]. Above hc2 the canted antiferromagnetic phase 



takes place, in wich the spins rotate according to the in- 
creasing field so that at /ic ~ 13 T they are aligned and 
the uniform, partially magnetized phase appears. For 
/i||[100], however, the canted ferrimagnetic phase is ab- 
sent, and at finite field only the canted antiferromag- 
netic and partially magnetized phases can be observed. 
Changing the polarization-polarization coupling to ferro 
type, the situation would be reversed: in ft,] [[110] there 
would be no sign of the canted ferrimagnet, while it would 
appear for ft.||[100]. Using the model that qualitatively 
recovers the [110]-field measurements, we gave a predic- 
tion to the field dependence of polarization when h\\ [100] . 
Applying finite temperature mean field theory we deter- 
mined the polarization for both orientation of magnetic 
field at various temperatures. The mean field results 
capture qualitatively well the temperature dependence 
of polarization and magnetisation: while the former is 
very sensitive to the temperature change, the latter is 
almost unaffected by it. Based on our calculations, we 
believe that relevant information regarding the low field 
phase could be obtained from torque measurements. Fur- 
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FIG. 8: Staggered polarization and magnetization results at 
finite temperature. We predict that in this rotated field set- 
ting the canted ferrimagnetic phase is missing, and the lower 
polarization curve is selected when > 0. Similarly to the 
h\ I [110] case, the polarization depends strongly on the temper- 
ature, however the magnetization hardly changes. (J = Jz = 
1.885 K, A = 15.08 K, = -0.02 K, = 0.01 K and 
g = 2.2) 
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